Coulomb effects on the formation of proton halo nuclei 
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The exotic structures in the 2sx/2 states of five pairs of mirror nuclei 17 0- 17 F, 26 Na- 26 P, 27 Mg- 27 P, 
28 A1- 28 P and 29 Si- 29 P are investigated with the relativistic mean-field (RMF) theory and the single- 
particle model (SPM) to explore the role of the Coulomb effects on the proton halo formation. The 
present RMF calculations show that the exotic structure of the valence proton is more obvious than 
that of the valence neutron of its mirror nucleus, the difference of exotic size between each mirror 
nuclei becomes smaller with the increase of mass number A of the mirror nuclei and the ratios of the 
valence proton and valence neutron root-mean-square (RMS) radius to the matter radius in each 
pair of mirror nuclei all decrease linearly with the increase of A. In order to interpret these results, 
we analyze two opposite effects of Coulomb interaction on the exotic structure formation with SPM 
and find that the contribution of the energy level shift is more important than that of the Coulomb 
barrier for light nuclei. However, the hindrance of the Coulomb barrier becomes more obvious with 
the increase of A. When A is larger than 34, Coulomb effects on the exotic structure formation will 
almost become zero because its two effects counteract with each other. 

PACS numbers: 21.10.Gv, 21.60.Jz, 21.10.Dr 



I. INTRODUCTION 

Nuclear halo is a kind of exotic structures in which the nuclear matter distribution extends to large radii because 
of its weakly bound character. Halo nuclei have been extensively investigated both experimentally and theoretically 
for decades 0, S i 1, i, i 0, S, I, E, E E E El E E E E E, E HI Hi- However, up to date most 
of the halo nuclei confirmed are the neutron halo. The proton halos observed are rather scarce. Since 8 B and 17 Ne 
are predicted as proton halos [IE HI, HE HI], the first excited state of 17 fE,H3 and proton- rich isotopes P and S 
[3 EE HE HE HE] are predicted as the proton halos and some of them are probed experimentally. 

There is a popular opinion that it is more difficult to form proton halo because the Coulomb barrier hinders the 
proton penetrating into the out region of nucleus. However, this may be not true actually in some cases. As will 
be seen below, for lighter nuclei the proton halo is easier to occur as compared to neutron halo due to the Coulomb 
interaction. The Coulomb interaction has two effects on the formation of nuclear halo. One of these effects is that it 
makes the energy level shift closer to the Fermi level thus facilitating the penetration of the valence proton beyond the 
range of nuclear force. On the other hand, the Coulomb barrier of the proton hinders the formation of halo structure. 
In the case of lighter nuclei, the former effect is more important than the later one. 

Because nuclear force is nearly charge-independent, the structure difference between mirror nuclei should mainly 
come form the effects of the Coulomb interaction. Therefore, the mirror nuclei in the neutron-rich side would act 
as a useful reference system to explore the role of the Coulomb interaction on the exotic structure formation. For 
this purpose, in this work we investigate the structures of five pair of mirror nuclei 17 0- 17 F, 26 Na- 26 P, 27 Mg- 27 P, 
28 A1- 28 P and 29 Si- 29 P which are all in the 2si/ 2 state in the framework of RMF. Combining with the RMF results of 
the 2si/2 states of 15 0- 15 N and 21 Ne- 21 Na [2I(, we find that for the seven pairs of mirror nuclei, the exotic structure 
of proton-rich nucleus is more obvious than that of its mirror nucleus, the difference of exotic size between the mirror 
nuclei becomes smaller with the increase of mass number A of the nucleus, both the ratios of the valence proton and 
the valence neutron root-mean-square (RMS) radius to the matter RMS radius in each pair of mirror nuclei decrease 
linearly with the increase of mass number A. In addition, in order to interpret the results above, we calculate in detail 
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two opposite effects of energy level shifted and the Coulomb barrier hindrance on the exotic structure formation by 
means of the single particle model (SPM) and find that the contribution of the energy level shift is more important 
than that of Coulomb barrier when A is small. However, the hindrance of Coulomb barrier becomes more obvious 
with the increase of A. When A is larger than 39, Coulomb effects on the exotic structure formation will almost 
become zero because its two effects counteract with each other. 



II. THE RMF RESULTS AND DISCUSSIONS 



As the relativistic mean-field (RMF) theory is a standard method for describing prop erties of the s phe rical nuclei 
and some deformed nuclei and its details can been found elsewhere such as Refs. [3l|, |32|, [33J, 12J, |35|, [36|, [37], [38|, [39j, S3, 
EH, El, El, S3, EE EE E3] , nere we om y describe the outline of the theory. For a system with the interacting nucleons, 
a, u>, and p mesons and photons, the Lagrangian density is written as 

C = ^[i^d^ fix g a a - g^w^ - 9 p Yt ■ 



■■^-A^ + \d»ad p o \ 



\r ■ + ■ P, - \A^A, V , (1) 



where 



= ^"-9V, (2) 
A itu = 9 » A v -Q V A», (3) 
r = d»(T - d v p» - 2g p (p» x ff), (4) 

The nucleon field and the rest mass are denoted as ip and m, respectively. The meson fields and their masses are 
denoted by <r, lo, p and m a , m u , m p , respectively. The photon field A p produces the electromagnetic interaction, and 
e is its coupling constant. g a , g u , g p are the coupling constants between the mesons and nucleons respectively, 32 and 
173 are the non- linear coupling constants of the a meson. T3 is the third component of the isospin Pauli matrices, i.e., 
r 3 \ n ) = \ n ) and r 3 \p) — ~~ \p)- Using Euler-Lagrang equation, a set of coupled equations for nucleons, mesons and 
photons can be obtained from the Lagrangian density function, which can been solved self-consistently by iteration 
under the mean-field approximation. After the final solutions are obtained, some quantities we need, such as the 
binding energy, single-particle levels, root-mean-square (RMS) radii of neutron and proton density distributions, and 
so on, can be calculated from the wave functions. In the effective lagrangian density function of the relativistic 
mean- field m is the mean value of the rest masses of proton and neutron, e satisfies e 2 /47r = 1/137 and rn CT , m w , 
m p, 9a, g^j, gp, 52, ff3 are free parameters. There are several well-tested nonlinear RMF parameter sets NL1, NL2, 
NL3 and NL-SH which were obtained by fit ting the experimental observables, such as the binding energies and radii 
of the nuclei HH, M, M, H S EE [H H HI S3, S3, E3, El, E3, El, EE E3- For usual RMF calculation there are 
two methods using spherical coordinate system and cylindrical coordinate system, which are used in calculations of 
spherical nuclei and reformed nuclei, respectively. 

In this work five pairs of mirror nuclei 17 F- 17 0, 26 P- 26 Na, 27 P- 27 Mg, 28 P- 28 A1 and 29 P- 29 Si and their core nuclei 
16 0, 25 Na, 25 Si, 26 Mg, 26 Si, 27 A1, 27 Si, 28 Si are calculated by RMF. In the calculations the parameter set NL1 (Mill 
is chosen for 16 0, 17 F, and 17 0, and NL3 [U EE] are chosen for other nuclei, and the Pauli blocking effects are 
considered. In addition, all nuclei are considered as spherical, which means the deformation effect is omitted to 
simplify the calculations. 



The RMF calculation results for the 2si/ 2 states of five pairs of mirror nuclei and their relative core nuclei are listed 
in Table |TJ £>exp and -B^he are ^ ne experimental and the calculated binding energies, respectively. R n , R p , R m and 
Rln denote the calculated RMS radii of neutron, proton, matter and the last valence nucleon density distributions, 
respectively, Eln is the single-particle energy of the last nucleon, Rln / Rm is the ratio of the valence nucleon RMS 
radius to the matter radius, and \R n — R p \ represents the difference between the proton and the neutron RMS radii. 
In Table [H B exp are taken from Ref.[48], in which B exp of the nucleus 26 P is the estimated value (denoted as #) by 
Audi and Wapstra because its experimental binding energy is unknown. The density distributions of neutron, proton, 
matter and the valence nucleus in the mirror nuclei 17 F- 17 which are calculated by the RMF, are shown in Fig.l. 
Here the normalized density distributions p{r), which satisfies / +O ° p{r)r 2 dr = 1, is shown in order to compare the 
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relative distributions of protons and neutrons. Based on Table [J and Fig. 1, the following discussions are made and 
some important conclusions are obtained. 

First, for the double magic nucleus 16 and its neighbor nuclei 17 and 17 F, the difference between the theoretical 
binding energy B t he and the corresponding experimental value B exp is very small, and B t he is at most 0.5% off. For 
other four pairs of mirror nuclei and their core nuclei, there are good agreement between their B t h e and B expi in 
which for four pairs of mirror nuclei, B t he is at most 2% off, and for other core nuclei B t he is at most 2.5% off. These 
results show that it is reasonable and reliable for us to use the spherical RMF theory to describes the properties and 
structures of the nuclei considered. 

Second, it can be seen from Fig. 1 that there are long tails in the density distribution of the valence proton of 17 F 
and the valence neutron of 17 0, compared with their core nucleus. The similar distributions also exist in other four 
pairs of mirror nuclei. These show that all of the five pairs of the mirror nuclei have exotic structures of halo or skin. 

Third, there are obvious differences in exotic size of each pair of mirror nuclei. For example, the binding energy eln 
of the valence proton is much lower than that of the valence neutron in its mirror nucleus, and the RMS radium Rln 
and \R p — R n \ of the valence proton is larger than that of the valence neutron in its mirror nucleus. In addition, it 
can also be seen from Fig.l that the tails in the density distribution of the valence proton of 17 F are longer than that 
of the valence neutrons of its mirror nuclei 17 0. The similar results occur in other four pairs of mirror nuclei. These 
show the exotic structure of valence proton is more obvious than that of the valence neutron of its mirror nucleus. 

Fourth, it can be seen from the change of the values Rln /Rm of the five pairs of mirror nuclei in Table Q] that the 
values Rln/Riti of the valence protons and the valence neutrons in each pair of the mirror nuclei all decrease with 
increase of mass number A. Also, the difference of exotic size between each pair of mirror nuclei becomes smaller 
with the increase of mass number A of the nucleus. In order to compare the size of the exotic structure between 
mirror nuclei more clearly, we plot Rln /Rm as a function of the mass number A in Fig. 2. The RMF results in 
the 2si/2 state of the other two pairs of mirror nuclei 15 0- 15 N and 21 Ne- 21 Na [2l[ are also shown in this figure for 
comparison. It can be seen from Fig. 2 that the values of 15 N, 17 F, 21 Na, 26 P, 27 P, 28 P and 29 P are consistently 
larger than those of their respective mirror partners, and the values RLN/Rm of the valence-proton nuclei and the 
valence-neutron nuclei are almost in two different lines and decrease gradually with the increase of mass number A. 
The two linear functions obtained with the least square fit are (Rln / Rm)p — 2.16176 — 0.02824A for the valence 
proton nuclei (valence proton line) and (Rln / Rm)n = 1-95596 — 0.02302/1 for the valence neutron nuclei (valence 
neutron line), respectively. The difference between the mirror nuclei decreases with the increase of mass number A, so 
as two lines intersect at about A=39. In addition, it is clear that the difference between the valence proton line and 
the valence neutron line also decreases linearly with increase of mass number A. This difference line (RMF Coulomb 
line) corresponds to the contribution of the pure Coulomb effect on the valence proton dispersion. 
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FIG. 1: The density distributions of neutron, proton, matter, and the last nucleon in 17 0- 17 F in the first excited state. Solid, 
dashed, dotted, and dash-dotted curves are the density distributions of protons, neutrons, matter and halo proton, respectively. 
The density distributions of their core nucleus 16 O is also drawn for comparing. 



The above results can be explained qualitatively by the RMF calculation results. Because the nuclear force is 
approximately charge-independent and there are the same nuclear force interaction in each pair of mirror nuclei, the 
difference between the structures of mirror nuclei should only come from the Coulomb interaction. The Coulomb 
effects on valence proton distribution include two opposite actions. One action is that it makes the energy level shift 
closer to the Fermi level, which means the valence proton is easier to extend far from the nuclear potential. The 
other action is that the Coulomb barrier hinders the diffusion of the valence proton and formation of exotic structure. 
We plot the single nucleon energy levels of 2si/ 2 states of five pairs of the mirror nuclei in Fig. 3 and the nuclear, 
the Coulomb and the total potentials for 17 0- 17 F and 29 Si- 29 P in Fig. 4. The potentials for the other three pairs of 
mirror nuclei are similar to that of 17 0- 17 F and 29 Si- 29 P. It can be seen from Fig. 3 that the single proton energy levels 
in 17 F, 21 Na, 26 P, 27 P, 28 P and 29 P are higher than the corresponding single neutron energy levels in their mirror 
nucleus 17 0, 26 Na, 27 Mg, 28 Al and 29 Si, which is induced by the Coulomb interaction. It can also be seen from Fig. 3 
that the energy level shifts in the 2s 1 / 2 states between mirror nuclei become larger with the increase of mass number 
A except 17 and 17 F. On the other hand, the calculation results shown that the hights of Coulomb barrier increase 
with the increase of mass number A. The Coulomb barriers of 17 F and 29 P are shown in Fig. 5. The Coulomb barrier 
of 26 P, 27 P, and 28 P lie between those of 17 F and 29 P (not shown in the figure). The two opposite effects of Coulomb 
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FIG. 2: The ratios of the last nucleon RMS radius to the matter one as a function of mass number in mirror nuclei. The solid 
stars stand for the ratios in the proton-rich nuclei and the solid squares stand for that in the neutron-rich nuclei. The results 
for 15 0- 15 N and 21 Ne- 21 Na are taken from Ref. [2l(]. The lines are the least square fit for the data. 
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FIG. 3: The single neutron energy levels of 2si/2 states for 17 0, 26 Na, 27 Mg, 28 Al, 29 Si and the single proton energy levels 
2s 1/2 states for 17 F, 26 P, 27 P, 28 P, 29 P. 



interaction determine the contribution of Coulomb potential to the valence proton dispersion together. Although the 
RMF theory cannot separate quantitatively the two opposite effects of Coulomb interaction to valence proton, it can 
be inferred from Fig. 4 and 5 that when mass number A is small, the contribution of the energy level shift is more 
important than that of the Coulomb barrier, but the hindrance of the Coulomb barrier becomes more obvious with 
the increase of A. When A is about 39, Coulomb effects on the exotic structure formation will almost become zero 
because its two effects counteract with each other. This qualitative conclusion will be proven quantitatively in the 
following single-particle model (SPM) analysis. 



6 




R(fm) 

FIG. 4: The variations of the means-field potentials with the radial coordinates for 17 0- 17 F and 29 Si- 29 P. Short-dashed, dashed 
and solid curves are the Coulomb potential, nuclear potential and total potential. 



TABLE II: The contributions of the energy level shifted and the Coulomb barrier to exotic formation calculated with the 
single-particle model. 

17 17 F 21 Ne 21 Na 26 Na 26 P 27 Mg 27 P 28 Al 28 P 29 Si 29 P 
Rln (fm) 4.27 5.43 4.40 5.37 4.19 4.87 4.10 4.60 3.98 4.33 3.93 4.22 
Ri (fm) 6.68 9.49 5.10 2.23 1.24 0.98 

R d (fm) -5.52 -8.52 -4.42 -1.73 -0.89 -0.39 
Rc (fm) L16 097 068 O50 0.35 0.29 



III. THE SPM RESULTS AND DISCUSSIONS 



In the single-particle model, the nucleus is assumed to be composed of the nuclear core and the valence nucleon 
outside the core. The normalized single-particle radial wave function in the (nlj) bound state 4> n ij(r) can be obtained 
by solving the Schrodingcr equation. The potential is chosen as 



V(r) = V N (r) + V c (r), 



(5) 
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FIG. 5: The Coulomb barriers of 17 F and 29 P. The solid and dashed curves are for 17 F and 29 P, respectively. 



where the nuclear potential Vpf(r) is chosen as Woods-Saxon potential 



and Vc(r) is Coulomb potential. Here the potential depth Vq of Woods-Saxon potential is adjusted to reproduce the 
valence nucleon separation energy, ro and ao are the radius and diffuseness parameter, respectively. In this work, ro 
and ao are chosen as the normal values 1.25 fm and 0.65 fm, respectively. The root-mean-square (rms) radius of the 
distribution of the last nucleon can be obtained from the single-particle radial wave function n y (r) by 



The numerical results with the SPM are listed in Table HH where Ri denotes the increase of Rln induced by the 
energy level shifted and Rd denotes the decrease of Rln induced by the Coulomb barrier. R c denotes the contribution 
of total Coulomb effect on the valence proton dispersion and R c = Ri + Rd- We select 17 0- 17 F as an example to 
explain the details of the calculations. At first, we obtain the RMS radii of the valence nucleon with the single-particle 
model (4.27 fm for 17 and 5.43 fm for 17 F). Next, we assume that the binding energy of 17 is the same as 17 F and 
calculate the valence neutron RMS radii for 17 under this condition. The result is 10.95 fm, which is much larger 
than the actual value of 17 0. This implies that the energy level shifted makes the RMS radius increase 6.68 fm if 
the effects of Coulomb barrier switch off. However, the Coulomb barrier hinders the exotic formation in reality. The 
RMS radius of 17 F is less than the result obtained by neglecting the effect of the Coulomb barrier. Therefore, the 
Coulomb barrier makes RMS radii decrease 5.52 fm. In Table HT1 "-" of Rd values shows decrease of the valence proton 
dispersion. It shows that two effects of the Coulomb interaction are opposite on the formation of exotic structures of 
valence proton, and the contributions of the energy level shifted are more important than that of the Coulomb barrier. 
The same calculations are applied to the other pairs of mirror nuclei. In order to show these rules more clearly, the 
variations of Ri, Rd and R c with the mass numbers are given in Fig. 6. It is seen from Table II and Fig. 6 that the 
two effects of the Coulomb interaction on the valence proton dispersion are opposite, and when mass number A is 
smaller the contribution of the energy level shifted is more important than that of the Coulomb barrier hindrance. 
However, the hindrance of the Coulomb barrier becomes more obvious with the increase of A, and total effect of the 
Coulomb interaction on the valence proton dispersion decreases linearly with increase of the mass number A of the 
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FIG. 6: The contributions of the energy level shift and the Coulomb barrier to exotic formation calculated with the single- 
particle approach. The solid circles denote the increase of Rln induced by the energy level shifted, the solid squares denote 
the decrease of Rln induced by the Coulomb barrier and the solid stars denote the total Coulomb effect on Rln- 



nucleus, and when A is larger than about 34, Coulomb effects on the exotic structure formation will almost become 
zero because its two effects counteract with each other. These conclusions are in good agreement with those obtained 
with RMF although there are some differences in their values, which maybe come from the difference between the 
two theoretical models RMF and SPM. 



IV. SUMMARY 



In summary, we have investigated the exotic structures in 2si/ 2 state of five pairs of mirror nuclei 17 0- 17 F, 26 Na- 26 P, 
27 Mg- 27 P, 28A1- 28 P and 29 Si- 29 P with RMF and SPM in order to explore the role of the Coulomb interaction on the 
proton halo formation. By analyzing the RMF results, we find that the exotic structure of valence proton is more 
obvious than that of the valence neutron of its mirror nucleus, the Rln I Rm values of the valence proton and the 
valence neutron of each pair of mirror nuclei decrease linearly with the increase of mass number A of the mirror nuclei, 
and the difference between the values of each pair of mirror nuclei becomes smaller linearly with the increase of mass 
number A. When A is about 39, the difference closes zero. By analyzing quantitatively two opposite effects of the 
Coulomb interaction on the valence proton dispersion in some mirror nuclei with SPM, we find that the contributions 
of the energy level shift are more important than that of the Coulomb barrier, and the total Coulomb effect on the 
valence proton dispersion becomes smaller linearly with increase of the mass number A, and it closes zero when A is 
larger than 34, which means two effects of Coulomb interaction on valence proton counteract with each other. These 
results basically agree with the conclusions obtained by RMF. 



9 



ACKNOWLEDGMENTS 

This work is supported by the National Natural Science Foundation of China under Grant Nos 10604008 and 
10435020, A foundation for the Author of National Excellent Doctoral Dissertation of China, and Beijing Education 
Committee Under Grant No. XK100270454. 



[1] I. Tanihata, H. Hamagaki, O. Hashimoto, Y. Shida et al, Phys. Rev. Lett. 55, 2676 (1985). 
[2] W. Mittig, J. M. Chouvel, W. L. Zhan, L. Bianchi et al, Phys. Rev. Lett. 59, 1889 (1987). 
[3] P. G. Hansen and B. Jonson, Europhys. Lett. 4, 409 (1987). 
[4] Y. Suzuki and K. Ikeda, Phys. Rev. C 38, 410 (1988). 
[5] M. G. Saint-Laurent et al, Z. Phys. A 332, 457 (1989) 
[6] E. Liatard et al, Europhys. Lett. 13, 495 (1990). 

[7] L. Johannsen, A. S. Jensen, and P. G. Hansen, Phys. Lett. B 244, 357 (1990). 

[8] G. F. Bertsch and E. Esbensen, Ann. Phys. (N. Y) 209, 327 (1991). 

[9] A. C. Hayes, Phys. Lett. B 254, 15 (1991). 
[10] M. V. Zhukov, B. V. Danilin, D. V. Fedorov et al, Phys. Rep. 231, 151 (1993). 
[11] I. J. Thompson, J. S. Al-khalili, J. A. Tostevin, and J. M. Bang, Phys. Rev. C 47, 1364 (1993). 
[12] Z. Z. Ren, B. Q. Chen, Z. Y. Ma, and G. O. Xu, Phys. Rev. C 53, R572 (1996). 
[13] Z. Z. Ren, W. Mitting, and F. Sarazin, Nucl. Phys. A 652, 250 (1999). 
[14] Z. H. Liu, C. J. Lin, H. Q. Zhang et al, Phys. Rev. C 64, 034312 (2001). 
[15] R. Morlock et al, Phys. Rev. Lett. 79, 3827(1997). 
[16] V. Guimaraes, S. Kubono et al, Nucl. Phys. A 588, 161c (1995). 
[17] L. Axelsson, M. J. G. Borge et al, Phys. Rev. C 54, R1511 (1996). 
[18] S. Aoyama, Phys. Rev. C 62, 034305 (2000). 

[19] N. Auerbach and N. Vinh Mau, Phys. Rev. C 63, 017301 (2000). 

[20] H. Y. Zhang, W. Q. Shen, Z. Z. Ren, Y. G. Ma et ai, Eur. Phys. J. A 23, 11 (2005). 

[21] J. G. Chen, X. Z. Cai, W. Q. Shen, Y. G. Ma et al, Eur. Phys. J. A 23, 11 (2005). 

[22] Y. J. Liang, X. H. Li, F. G. Deng et al, Chin. Phys. Lett 23, 1142 (2006). 

[23] K. Riisager and A. S. Jensen, Phys. Lett. B 315, 6 (1993). 

[24] J. S. Al-khalili, J. A. Tostevin, and I. J. Thompson, Phys. Rev. C 54, 1843 (1996). 
[25] H. Esbensen, Phys. Rev. C 53, 2007 (1996). 
[26] B. J. Cole, Phys. Rev. C 58, 1866 (1998). 

[27] Z. Z. Ren, A. Faessler, and A. Bobyk, Phys. Rev. C 57, 2752 (1998). 
[28] B. A. Brown and P. G. Hansen, Phys. Lett. B 381, 391 (1996). 
[29] B. Q. Chen, Z. Y. Ma et al, J. Phys. G: Nucl. Phys. 24, 97 (1998) 
[30] A. Navin et al, Phys. Rev. Lett. 81, 5089 (1998). 
[31] J. D. Walecka, Ann. Phys. 83, 491 (1974) . 

[32] C. J. Horowitz and B. D. Serot, Nucl. Phys. A 368, 503 (1981) . 
[33] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 (1986) . 

[34] P. G. Reinhard, M. Rufa, J. A. Maruhn, W. Greiner, and J. Friedrich, Z. Phys. 323, 13 (1986). 

[35] M. Rufa, P. G. Reinhard, J. A. Maruhn, W. Greiner, and M. R. Strayer, Phys. Rev. C 38, 390 (1988). 

[36] P. G. Reinhard, Rep. Prog. Phys. 52 439 (1989). 

[37] I. Tanihata, D. Hirata, T. Kobayashi, S. Shimoura, K. Sugimoto, and H. Toki, Phys. Lett. B 289, 261 (1992). 
[38] D. Hirata, H. Toki, 1. Tanihata, and P. Ring, Phys. Lett. B 314, 168 (1993). 
[39] L. S. Warrier and Y. K. Gambhir, Phys. Rev. C 49, 871 (1994). 
[40] P. Ring, Prog. Part. Nucl. Phys. 37, 193 (1996). 

[41] G. A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C 55, 540 (1997). 
[42] B. D. Serot and J. D. Walecka, Int. J. Mod. Phys. E 6, 515 (1997). 

[43] J. Meng, K. Sugawara-Tanabe, S. Yamaji, P. Ring, and A. Arima, Phys. Rev. C 58, R628 (1998). 
[44] J. Meng, K. Sugawara-Tanabe, S. Yamaji, and A. Arima, Phys. Rev. C 59, 154 (1999). 

[45] Z. Patyk, A. Baran, J. F. Berger, J. Decharge, J. Dobaczewski, P. Ring, and A. Sobiczewski, Phys. Rev. C 59 , 704 (1999). 
[46] M. Bender, P. H. Heenen, and P. G. Reinhard, Rev. Modern Phys. 75, 121 (2003). 
[47] Y. S. Li and G. L. Long, Commun. Theor. Phys. 41, 429 (2004). 
[48] G. Audi and A.H. Wapstra, Nucl. Phys. A 565, 1 (1993). 



